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SIMULATION  OF  UNDERKATE?  DIVERSITY  ARRAY 

In  earlier  reports  [1,2]  a  three-cinons ion.il  underwater  array  which  simul¬ 
taneously  forr.s  and  combines  multiple  beams  in  elevation,  vas  described  and! 
analysed.  •  The  objective  of  tbo  system  is  to  approve.'.,  a  condition  in  vhicL  the 
multiple  ray  arrivals  from  a  distant  source  r. r o  separately  received  and  coherently 
combined-  Because  the  analytical  forms  giving  t.hc'  final  output  signal -tc-noise 
ratio  (AVI;)  are  too  involved  for  direct  ccnputotion^  s  L.. elation  experiments  mere 
carried  out.  Results  of  thin  work  are  presented  fee l aft,  for  ttuV  system  described 
in  [1  ,  section  3,  {»p.  26-3H  v;hi  Hi  js  bas*-*d  on  niirdr.al  ra  t  vo  'ennj  i  w  inf*  (MKC)  of 
the  r.'iill Ip  1  i :  rays  for:r»e»i.  At  the  same  lime  certain  other  cases-  vrere  simulated  ^ 
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in  order  to  see  what  magnitude  of  improvement  is  obtained  using  this  method  over 
simpler  ones.  Given  below  arc  results  for  (1)  a  three-dimensional  array  with 
multiple  beams  (branches)  using  only  coherent  phasing  of  branches  (known  as 
equal  gain  combining  (EGC) ,  (2)  a  three-dimensional  array  with  multiple  beams 
using  selection  of  the  maximum  amplitude  branch  (known  as  selection  combining 
(SC)),  (3)  a  three-dimensional  array  with  a  single  fixed  focus  beam*  and  (A)  a 
two-dimensional  arrsy.^The  last  one  represents  the  original  concept  explored 
earlier  and  reported  a/n  "{3,4,51  . 

The  simulation  results  can  be  brief l}r  summarized  as  follows.  For  the  con¬ 
ditions  chosen,  an  improvement  of  approximately  A  to  1  is  obtained  using  the 
maximal  ratio  combining  technique  with  the  three-dimensional  array  over  anything 
else  that  vras  done.  An  interesting,  though  expected  result,  was  that  the  var¬ 
iance  relative  to  the  mean  or  the  output  SNR  using  this  best  technique  was  much 
less  than  for  the  two-dimensional  array.  The  computational  model  assumed  ten 
independent  ray  arrivals  randomly  spread  over  ±10°  in  the  vertical.  In  the  two- 
dimensional  array  those  combine  noncoherent lv  resulting  in  a  nearly  Rayleigh 
fluctuation  of  amplitude.  In  the  three-dimensional  array  the  ten  rays  are  to 
a  large  extent  resolved  in  the  separate  branches.  The  diversity  selection  or 
coherent  combining  then  results  in  a  substantially  smaller  fluctuation. 

We  now  describe  the  physical  and  statistical  arrangements  assumed.  Fig¬ 
ure  2.1  suggests  the  deployment  of  elements  in  a  three-dimensional  space.  Element 

positions  (X  ,Y  ,Z  } ,  n  -  1,2  N  were  assumed  independent  random  vectors,  the 

n  n  a 

number  of  elements  N  being  31  for  this  computation.  The  horizontal  coordinates 
<V  Y  )  were  assumed  independent  normally  distributed  random  variables  with  zero 
mean  and  standard  deviation  of  50  wavelengths.  This  value  of  standard  deviation 
impl  ies  that  at  100  Kz  where  the  wavelength  is  about  15  meters,  about  68%  of  the 
array  elements  will  be  concentrated  in  a  tango  of  t  750  meters  around  the  center 
of  the  array.  Because  the  array  main  beam  wan  focused  to  look  in  the  Y-Z  plane 
only,  the  random  variables  were  not  involved  in  the  computation.  The  vertical 
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FIGURE  2.1  PHYSICAL  ARRAN G EM ENT  ASSUMED  IN  SJHULATION. 


coordinate  Z  was  assumed  unifonr.lv  distributed  with  mean  zero  and  range  of 

ri 

50  wavelengths,  also.  This  implies  a  depth  range  of  *  375  meters  around  the 
array  center. 

Ray  arrivals  were  assumed  to  be  in  the  Y-Z  plane  corresponding  to  the  azi¬ 
muthal  angle  of  focus  of  the  array,  but  the  rays  were  assumed  dispersed  in 
vertical  angle.  Ten  rays  were  assumed  arriving,  all  of  equal  magnitude  and 

independent  random  electrical  phases  :t  ,  m  -  1,2  ...  10,  each  uniformly  d  is- 

m 

tributed  in  (0.  2~i).  The  arrival  angles  0  ,  measured  from  the  vertical  were 
assumed  independent  and  uniformly  distributed  over  (80°,  100°),  or  i  10°  rela¬ 
tive  to  the  horizontal. 

Ten  sets  of  random  pairs  of  numbers  (y  ,z  )  n  ~  1»2  ...  31  to  represent 

ten  possible  random  element  positions  were  chosen.  For  each  of  these  positions, 

five  sets  of  angle  pairs  (3  ,ct  )  ,  m  r*  1,2  ...  10  were  randomly  selected. 

mm 

Assuming  the  combining  scheme  of  [1,  section  3,  Figuie  3.2],  the  output  of 
ca~h  branch  prior  to  weighting  and  combining  is  given  by 
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where  8  )  are  the  azimuth  and  elevation  angles  of  the  mth  arriving  ray,  and 

(B  ,  a  )  are  the  magnitude  and  electrical  phase  at  the  center  of  coordinates  of 
tu  rn 

that  ray.  is  the  injected  electrical  angle  to  focus  the  array  to  azimuth  and 

elevation  angles  (4>g,  ©si)  as  shown  in  [1]  Figure  3.2,  page  31,  and  is  given  by 

4>  .  K  -k[x  sin0  .  cos$  +y  sinB  *sind>  -fz  cosO  A  (2) 

Tni  n  si  s  Jn  si  n  si  v  7 


The  subscript  i  identifies  the  beam  or  branch  number  and  I  is  the  number  of  branches 

used;  here  we  will  use  I  =  17  branches  for  reasons  to  be  stated  below.  (1)  and  (2) 

are  written  to  include  rays  arriving  from  any  azimuth  and  elevation  (a  ,  ©  )  and 

focused  at  any  azimuth  and  elevation  ($ g>  0 g^) .  We  now  specialize  these  expressions 

to  the  case  of  the  array  focused  in  the  Y-Z  plane  and  ray  arrivals  in  the  Y-Z  plane. 

Thus  <J  =  <f>  -  90°  and  (1)  and  (2)  reduce  to 
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The  final  signal  output  using  KRC  is  given  by 

s  -  \  iv  1 2  I  A2  (5) 

i«l  i~?.  1 

The  mean  square  value  of  the  noise  at  the  output  of  the  system  is  given  by 


N 

I 

n=l 


<  n 


I 

I 

i=l 


y 


-1 


A  .A. co. 
i  J 


(6) 


where  the  (A.,  . )  are  the  measured  signal  amplitude  and  phase  on  the  ith  branch 
x  1  2 

given  in  (1),  the  ^  are  given  by  (A),  and  <n^>  Is  the  equivalent  moan  square 
value  of  the  noise  input  genera t ed  by  the  nth  sensor.  In  t be  computation  <n^> 
was  set  equal  to  1.  Finally  the  computed  output  SNR  is  given  by 
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S2 

SNU  =  (7) 

N 

As  discussed  in  [2]  the  angles  0  . ,  if  set  at  values  separated  by 

s  x  n 

n  an  integer,  will  result  in  uncorrelated  noise  voltages  at  the  branch  output  pro¬ 
vided  the  elements  are  uniformly  distributed  in  depth  over  the  range  (-h,  h) . 

For  the  geometric  conditions  used  here  it  was  determined  that  1,15°  spacing  between 
the  0  ^  would  accomplish  this  when  h  =  25X.  We  have  accordingly  assumed  17  beams 
symmetrically  placed  around  6  =  90°  at  1,15°  intervals.  Separation  is  about  a 
beamwidth  in  this  case  and  the  total  coverage  in  vertical  angle  is  close  to  ±10° 
with  respect  to  the  horizontal. 

Results  of  the  simulation  using  the  maximal  ratio  combining  technique  are 
shown  in  Table  1.  Here  we  show  the  average  and  standard  deviation  of  the  SNR  given 
by  (7)  over  the  five  sets  of  paired  values  (0^,  a^)  m  -  1,  2,  ...  10,  for  each  of 
ten  sets  of  position  samples  (y^,  z^)  n  =  1,  2,  ...  31.  These  statistics  are  de¬ 
noted  <SNR>^a  and  a (SNR)...  Then  the  50  results  of  SNR  (five  sets  over  (0  ,  a  ) 

09  09  m  m 

times  the  ten  sets  over  (y^,  z^))  were  treated  as  a  sample  of  size  50  and  the 

overall  mean  and  standard  deviation  denoted  <SNR>  and  a  respectively  were  determined* 

•  • 

These  were  found  to  be  <SNR,>  =  480  and  a  -  182.  We  point  out  that  if  a  single  ray 
were  assumed  to  impinge  on  a  single  element  the  output  SNR  would  be  unity  and  the 
variance  would  be  zero. 

To  determine  the  effect  of  a  frequency  change  on  the  output  we  have  assumed 
two  different  situations.  In  the  first  we  assumed  the  angular  separation  between 
beams  held  at  1.15°:  that  is,  0  ^  was  held  fixed  at  90°  ±  n(1.15°),  n  -  0,  1,  2, 

...8.  The  frequency  was  then  changed  by  factoxs  1/2  and  2.  These  latter  changes 
were  accomplished  bv  simnly  changing  the  value5  of  y  and  z  used  in  the  previous 
calculation  by  the  reciprocal  of  these  same  factors.  The  vertical  beamwidths  be¬ 
come  narrower  at  the  higher  frequency  and  wider  at  the  lover  frequency  but  the 
angular  spacing  between  beams  remains  unchanged.  Thus  the  beams  are  not  optimally 
spaced  resulting  in  either  uncorrelated  noise  in  the  several  branches,  or  in  non- 
total  coverage  of  the  vertical  range  within  which  incoming  rays  are  expected.  Result! 
of  these  calcul at ions  are  also  shown  in  Table  1  revealing  a  decrease  in  the  over-  , 
all  average  SNR  as  one  might  expect.  These  results  essentially  shew  the  sensitivity 
of  the  scheme  to  incorrect  placement  of  the  vertical  beams.  As  x</e  see,  the  effects 
are  not  overly  serious,  the  mean  output  at  the  1/2  and  2  times  frequency  points 
being  within  about  80%  of  the  menu  at  the  design  frequency. 
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In  the  second  simulation  of  frequency  effect  the  array  focusing  was  held 
fixed  for  a  100  Hz  sinusoid  by  fixing  the  angles  The  applied  frequency  was 

then  altered  to  99,  99.5  and  100.5  Hz.  As  a  rule  of  thumb  the  bandwidth  of  an 
array  with  fixed  phase  shift  focusing  is  the  inverse  of  the  time  required  for 
the  wave  to  traverse  the  array.  In  this  case  it  would  imply  a  bandwidth  of  the 
order  of  1  Hz,  or  about  1%  of  the  center  frequency.  The  results  of  the  simula¬ 
tion  are  shown  in  Table  2.  The  overall  mean  SNR  is  observed  to  have  fallen  by 
about  3  dB  at  frequencies  100  ±0.5  Hz  from  what  it:  was  at  100  Hz,  thus  confirming 
the  rule  of  thumb  on  bandwidth.  The  overall  mem  SNR  at  99  Hz  has  fallen  further, 

the  level  appearing  to  be  about  what  one  gets  when  one  steers  the  azimuthal  focus 

& 

away  from  the  source,  illuminating  the  sidelobes. 

We  point  out  that  the  array  properties  observed  here  are  all  normalized  to 
wavelength  so  that  at  higher  frequencies,  with  the  actual  array  size  reduced  but 
with  array  size  in  wavelengths  held  constant,  the  bandv;idth  would  remain  at  about 
1%.  At  10  kHz  we  expect  a  100  Hz  bandwidth,  a  value  adequate  for  operating  a 
teletype  communication  link.  Furthermore,  at  this  center  frequency  the  array  hor¬ 
izontal  dimension  measured  between  to  points  is  15  meters,  a  dimension  one  might 
envision  for  an  array  suspended  from  a  surface  ship  or  deployed  around  a  submerged 
submarine.  The  array  could  therefore  he  useful  for  underwater  data  communication. 
Furthermore,  our  wark  was  based  on  focusing  by  fixed  phasing  of  elements  of  ;he 
array.  It  is  the  fixed  phasing  which  limits  the  array  bandwidth.  By  using  con¬ 
trolled  time  delay  networks  at  each  element  involved,  broader  bandv/idths  are  achiev¬ 
able  suggesting  the  possibility  of  higher  speed  data  communication,  or  lower  center 
frequency  with  higher  bandwidth. 

The  nultibcnn  three-dimensional  array  was  compared  to  four  other  arrangements 
as  follows: 

1.  The  maximal  ratio  combiner  weighting  circuits  which  multiply  each  branch  out- 
put  by  *=  A^e  J  i  where  V_.  is  given  by  (1)  are  replaced  by  constant*  amplitude 
phase  shifters  e  ^ri,  thus  (5)  and  (6)  become 

1 

s  »  l  A.  (8) 

i=l  1 


T)ic*  rsideloW  proven  ii-<;  , i r e  iHi'.citnsod  furHirr  below. 
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Here  only  phase  tracking  is  needed  but  the  mean  SNR  will  not  be  as  good  as  for 
maximal  ratio  combining.  The  term  "equal  gain  combining  (EGC) "  is  used  in  diver¬ 
sity  communication  for  this  arrangement. 

2.  The  maximal  ratio  combiner  was  replaced  by  a  "selection  combiner  (SC)H;  that 
is,  one  which  simply  selects  the  output  with  the  maximum  SNR.  This  technique  is 
also  a  standard  scheme  in  communication  diversity  systems.  Whereas  the  maximal 
ratio  combiner  gives  an  output  SNR  which  is  the  sum  of  branch  SNRs,  this  scheme  . 
produces  only  the  maximum  of  the  branch  SNRs.  There  is,  however,  no  need  for 
phase  tracking,  greatly  simplifying  the  processing.  It  is  possible  that 

in  applications  such  as  the  underwater  case  the  problem  is  as  much  variation 
of  vertical  arrival  angle  as  it  is  multipath.  In  effect  then,  the  array  would 
follow  the  variation  in  angle  of  arrival  of  the  maximum  amplitude  ray. 

3.  Output  was  taken  from  one  branch  of  the  multi-beam  array,  the  one  which  fo¬ 
cuses  horizontally  (r  .  =  90°).  The*  our;  ose  of  this  calculation  is  to  see  what 
effect  is  obtained  when  the  threc-dimas  j  onal  array  is  operated  in  its  simplest  mode. 

4.  The  three-dimensions]  array  was  reduced  tc  a  planar  horizontal  array.  Here 
we  are  returning  to  the  original  array  structuie  —  the  two  dimensional  arrey. 

Sumary  results  ‘rearing  out  expectations  r re  shown  in  Table  3.  The  maximal 
ratio  combining  sc  he  re  gives  an  overall  mean  SNR  of  at  least  6  dB  better  than  the 
other  arrangements  except  for  HOC  case  which  is  a  close  second.  Interestingly , 
the  ratio  of  c*  (SNR) /<SNR>  is  much  small  c r  in  the  diversity  nodes  than  in  the  two- 
dimensional  case.  This  too  is  expected.  The  diversity  inodes  arc  partly  effective 
in  resolving  the  multipath  and  avoiding  the  nou-cohercat  interference  of  the  multi- 
path  ccrf  o.it  ntfi;  In  the  two*  d  imension*'  l  case  all  rays  entering  the  relatively  wide 
vertical  array  boanwidth  are  combined  non-coherc ut  i y .  There  is,  therefore,  consid¬ 
erable  amplitude  '\irintion  d-g-nding  on  the  relative  phases  of  the  accepted  rays. 

In  the  fixed  bo..m  t  lire.- -d  r.  .•  nsionnl  rare  there  is  also  a  high  ratio  c  (SNR) /<  SNR> 
presiiiM.ibly  a  ro.uf?  of  the  narrow'  vertical  bear-width  which  may  or  may  not  see  ar¬ 
riving  acoustic  energy. 

An  important  property  ol  tin-  array  will  bo  its  response  to  sources  off  the 
azimuth  of  focus.  Tf  ve  wore  t<>  wing,  ic  swing  in*,  the  focus  away  from  a  source  which 
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MRC 

EGC 

SC 

90°  Sector 

2-Dimensional 

<SNR> 

480 

367 

105 

32 

117 

o (SNR) 

182 

156 

48 

30 

107 

o  =  50X  for  an  cases 

y 

h  =  25X  for  3-dimensional  arrays 

TABLE  3  COMPARISON  OF  SIMULATION  RESULTS 

neatly  fed  M  branches  independently  the  power  in  each  branch  would  on  average  drop  by 

a  factor  N,  the  number  of  elements.  But  the  final  output  with  MRC  or  EGC  being  the  < 

M  ! 

coherent  sum  of  the  M  weighted  branches,  would  only  be  reduced  by  a  factor  —  from 
the  main  beam  power.  A  situation  of  this  sort  was  simulated  to  check  this  surmise. 
The  result  obtained  using  the  MRC  system  v as  <SNR>  --  182,  the  average  being  over  the 
same  set  of  random  variables  as  before.  This  figure  is  somewhat  below  (tt/N)  (main- 
beam  <SNR>)  but  it  indicates  that  these  techniques  do  exact  a  price  in  sidelole 
response.  In  this  calculation  there  was  no  signal  on  the  mainbenm. 

Further  study  of  side! oho.  effects,  with  and  without  a  main  beam  signal.  nvcscr.t 
and  using  the  diife  cent  combining  schemes,  is  viewed  as  a  useful  next  step.  In 
addition,  no tliocs  bason  on  estimation  theoretic  principles  (e^g«,  maximum  likeli¬ 
hood  and  maximum  entropy  estimation)  should  be  considered  for  application  here. 

These  methods  inherently  maximize  on-target  sign*.]  response  relative  to  off-target 
signals.  Applied  to  the  separate  beams  as  found  here,  or  even  to  the  entire  array, 
superior  sidelobe  rejection  characteristics  can  be  expected. 

Fred  Haber 
Paul  Ych 
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